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ABSTRACT 

We present a detailed analysis of the Horizontal Branch of the Carina Dwarf Spheroidal Galaxy by means of synthetic modelling 
techniques, taking consistently into account the star formation history and metallicity evolution as determined from main sequence 
and red giant branch spectroscopic observations. We found that a range of integrated red giant branch mass loss values of 0.1-0.14 M© 
increasing with metallicity is able to reproduce the colour extension of the old Horizontal Branch. Nonetheless, leaving the mass loss 
as the only free parameter is not enough to match the detailed morphology of Carina Horizontal Branch. We then investigated the 
role played by the star formation history on the discrepancies between synthetic and observed Horizontal Branches. We derived a 
“toy” bursty star formation history that reproduces well observed horizontal branch star counts, and also matches qualitatively the 
red giant and the turn off regions. This bursty star formation history is made of a subset of age and [M/H] components of the star 
formation history based on turn off and red giants only, and entails four separate bursts of star formation of different strenghts, centred 
at 2, 5, 8.6 and 11.5 Gyr, respectively, with mean [M/H] decreasing from ~ -1.7 to ~ -2.2 when the age of the burst increases, 
and a Gaussian spread of cr 0.1 dex around these mean values. The comparison between the metallicity distribution function of our 
bursty star formation history and the one measured from the infrared CaT feature using a CaT-[Fe/H] calibration shows a qualitative 
agreement, once taken into account the range of [Ca/Fe] abundances measured in a sample of Carina stars, that causes a bias of the 
derived [Fe/H] distribution toward too low values. In conclusion, we have shown how the information contained within the horizontal 
branch of Carina (and dwarf galaxies in general) can be extracted and interpreted to refine the star formation history derived from red 
giants and turn off stars only. 

Key words, galaxies: dwarf - galaxies: evolution - galaxies: stellar content - Hertzsprung-Russell and C-M diagrams - stars: 
horizontal-branch 


1. Introduction 

Dwarf galaxies (DGs) play a major role in modern astrophysics 
as they are believed to be the building blocks of the process 
of galaxy formation. Therefore DGs constitute a sort of fossil 
record of the formation epoch of the cosmic structures, and the 
determination of their star formation history (SFH) is crucial to 
understand the mechanisms of galaxy formation and early evo¬ 
lution. 

Detailed SFHs of DGs can only be determined in the Local 
Group where they can be resolved into individual stars down 
to the oldest main sequence (MS). These determinations are 
usually based on the theoretical interpretation (via stellar evo¬ 
lution models and isochrones) of observed colour magnitude 
diagrams (CMDs) and, when available, spectroscopic heavy- 
element abundances, typically of the red giant branch (RGB) 
populations. The primary age indicators for these populations 
are found in the main sequence turn-off (TO) region of the CMD, 
which is located at faint magnitudes for the oldest populations, 
and is very sensitive to photometric errors. 

The horizontal branch (HB) is routinely neglected in SFH 
determinations of DGs, in spite of its brightness compared to 
the old TO, and the extreme sensitivity -in terms of colour 
and brightness distribution- to the mass and metallicity distri¬ 
bution of the parent stars. The reason is that the interpretation 
of the HB morphology in potentially simpler populations like 


Galactic globular clusters is problematic. Numerous studies of 
Galactic globular clusters (GCs) have shown that age and metal¬ 
licity alone can not account for the mean colour and extension 
of the HB (e.g., |Catelan|2009[ [Potter et al.|2010t|Gratton et al. 


|2010| ). A major difficulty is that stars arriving on the Zero Age 
HB (ZAHB) have lost mass during the previous RGB phase, and 
to date it is still impossible to predict from first principles the 
amount of mass lost by RGB stars. This issue is further com¬ 
plicated by the currently well established presence in individ¬ 
ual GCs of multiple populations of stars with enhanced helium 
abundances at fixed metallicity (e.g. jGratton et al.|2012| and ref¬ 
erences therein), which affect both the colour and magnitude of 
the HB. 


Despite these complications, and the presence of He- 
enhanced populations at hxed metallicity may well be just a fea¬ 
ture of GCs, HB stars contain a wealth of information to con¬ 
strain the star formation rate and metallicity evolution of DGs 
at earliest times. In jSalaris et al.j ( |2013| ) we presented the hrst 
detailed simulation of the HB of the Sculptor Dwarf Spheroidal 
(dSph) galaxy by means of synthetic modelling techniques, tak¬ 
ing into account the SFH and metallicity evolution determined 
from the MS and RGB spectroscopic observations. We found 
that the number count distribution along the observed HB could 
be reproduced with a simple ma ss loss law (that a grees very 
closely with the determinations by jQriglia et al.|20f4 for a sam- 
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Fig. 1. Carina dSph. Upper panel: Star Formation Rate as a 
function of time from the adopted SFH. Lower panel: Weighted 
mean metallicity as a function of time. The red dotted histograms 
show the Icr dispersion. 


pie of GCs), and that there is no excess of bright stars that require 
He-enhanced populations. 


The purpose of the present work is to investigate, by means 
of the same synthetic HB modelling as for Sculptor, the HB pop¬ 
ulations of another DG belonging to the Local Group, the Carina 
dSph galaxy. This galaxy has a SFH very different from that of 
Sculptor, showing multiple star formation episodes separated in 
time and in the amount of stellar mass involved. Numerous stud¬ 
ies have probed a broad range of properties, through deep pho- 


tometric investigations (e.g. Smecker-Hane et al.|1996t|Monelli| 

et al.||2003l |Bono et al.||20101); spectroscopic anal 

ysis, : 

)oth at 

medium (e.g.|Smecker-Hane et al. 1999[| 

Koch et al. 

'sO 

o 
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(N 

Jelmi 

et al.|20061) and high (e.g. Shetrone et al. 

120031 Koc 

1 et a . 

2008 

Fabrizio et al.||2012[ |Lemasle et al.||2012| |Venn et al. 

2012) 

resolution; variable star characterisation (e.g. Saha et al. 

1986 

Mateo et al. 1 1998 [ palT Ora et al.||2003 Coppola et al. 

2013) 

and SFH analysis (e.g. Pasetto et al.||2011[ [Small et al. 

2013 


|de Boer et al.||2014| ). These detailed works provide us with the 
information needed to make a meaningful comparison between 
synthetic and observed HB populations. In particular, our anal¬ 
ysis will provide strong additional constraints on the controver¬ 
sial issue of the galaxy metallicity distribution function (MDF) 
as determined spectroscopically and photometrically from RGB 
stars (see, e.g., |Bono et al.|2010[[VandenBerg et al.|2015l for dif¬ 
ferent conclusions about the consistency between spectroscopic 
and photometric MDFs). 


This work is structured as follows: we present the data set 
used for our investigations in §2; compare the synthetic and ob¬ 
served HBs in §3; discuss the impact of our results on the SFH 
in §4 while a summary is presented in §5. 


Fig. 2. Carina dSph. Comparison between the different fields of 
our da ta-set. The black dots represent stars of |de Boer et af 


( |2014| ) photometric catalog ue inside one tidal r adius. The blue 
rectangle is the field of the |Bono et ( [2010| ) catalogue. The 
red stars mark the position of RR Lyrae stars while t he blue 
squares are the stars with spectroscopic measurements ( [Helmi 


et al.|2006|). The zero point of the horizontal and vertical axis is 
set to RA = 6/z41m 36.6^ and Dec =-50° 57' 58" (J2000). 


2. Data 


For the computation of our synthetic HB we employed, as ref¬ 
erence SFH, that determined by |de Boer et al.|p014| ) (hereafter 
dB14). Among the several SFHs in literature, we chose the dB14 
solution because it is the only one that combines the photometric 
modelling of the CMD with spectroscopic information about the 
metallicity distribution. Due to the very high sensitivity of the 
HB morphology to metal content, this approach is preferable. 

For the sake of homogeneity we adopted the dB14 solution 
obtain ed with the BaSTI evolutionary tracks ( [Pietrinferni et al. 
2004), which are the same ones used in our synthetic HB calcu¬ 
lations. Carina photometry has been divided by dB14 into three 
concentrical annuli inside the tidal radius of the galaxy, plus a 
fourth field external to the tidal radius. For each field, an in¬ 
dependent SFH has been computed. The solution gives the star 
formation rate in a grid of age and [Fe/H] bins, with also an es¬ 
timate of [a/Fe] in each bin. The SFH used in our synthetic HB 
calculations refers to the sum of the annuli inside the tidal radius. 
Our Synthetic HB simulation with this SFH will be denoted as 
reference simulation. 

Figure shows star formation rate and mean stellar [Fe/H] 
as a function of time. There are at least two major epochs of star 
formation at old and intermediate ages, ranging between 9-14 
and 3-7 Gyr respectively, while the mean metallicity remains in 
the relatively narrow interval -2 < [Fe/H] < -1.6. 

We compare the synt hetic CMDs with the photometric data 
from |Bono et al.| ( [20101 ) (hereafter BIO), which include 4152 
CCD images acquired between December 1992 and January 
2005. Although the use of dB14 photometry would have granted 


2 




















































































































A. Savino et al.: Horizontal Branch stars and star formation history of Carina 



(B-V) 


our photometry, and replaced them with their intrinsic position 
on the CMD. We used the period distribution to add additional 
observational constraints to our HB simulations. 

A glance at the Carina (V,B- V) CMD (Fig. reveals the 
complexity of its stellar population. Two distinct MS Turn-Offs 
can be clearly seen, with distinct SGBs merging in a very nar¬ 
row RGB. The accepted scenario envisages Carina as a system 
that has undergone two or more major events of star forma- 
tion, separated by a period of several Gyr ( |Smecker-Hane et ah] 
1996t Monelli et al.|2003] ). In addition, a very young popula¬ 


tion is probably present, as suggested by the presence of a blue 
plume abo ve the MS and by th e detection of several Anomalous 
Cepheids ( [Monelli et al.|2003] ). The narrow RGB, coupled with 
the observed broad metallicity distribution ( Koch et 11] 2006[ 
Helmi et al.|2006| ), suggests a conspiracy among age, metal con¬ 


tent and alpha element abundances that leads to all stars having 
a very similar colour on the RGB. 

The burstiness of Carina SFH can be seen in the helium burn¬ 
ing loci as well, with an old extended Horizontal Branch that is 
clearly detached from a younger, more populated Red Clump 
(RC). The discretness of the Carina stellar populations is very 
helpful since it allows us to simplify our analysis and model each 
component separately. 


Fig. 3. CMD of the Carina dsph from the BIO photometric set 

a perfect match between the stellar population sampled and the 
SFH, we chose to employ a different data-set because of the sev¬ 
eral advantages it involves. 

First, the smaller photometric errors in the B10 CMDs, of the 
order of ~ 0.004 mag at the HB level, allow for a more robust 
comparison between observed and synthetic CMDs. In addition, 
the same field has been the subject of a deep variable star search 
(Coppola et al. submitted). Furthermore, as the BIO photometric 
catalogue merges observations taken at different epochs, differ¬ 
ent random phases in different frames, for a given variable star, 
tend to be averaged. The colour and magnitude of possible un¬ 
detected RR Lyrae in B10 will therefore be closer to the intrinsic 
mean values than it is for single epoch data-set. 

Finally, the photometric data-set from dB14 retains the con¬ 
tamination by foreground stars of the Milky Way, as it is statisti¬ 
cally taken into account at later stages during the SFH determi¬ 
nation. Since Carina is a fairly diffuse stellar system at relatively 
low Galactic latititude, the number of foreground stars is large 
and, because the HB is much less populated compared to other 
evolutionary stages, the contamination can be a serious problem. 
The BIO catalogue, instead, has been carefully cleaned of Milky 
Way stars and unresolved galaxies. 

Fig. shows the comparison between the dB14 field inside 
the tidal radius and the BIO field. 

In theory, a spatial gradient in the stellar population proper¬ 
ties may affect our comparison, due to the different areas sam¬ 
pled. Nonetheless, as the difference in the spatial coverage is 
a minor fraction of the total sampled area, and since the star 
counts are dominated by the central region of the galaxy, we ex¬ 
pect this effect to be quite small. It should be also noted that the 
results presented in this paper still hold when the synthetic CMD 
is compared with the dB14 photometric catalogue. 

For the variable star modelling, we used the catalogue from 
Coppola et al. (2015, submitted), which contains information 
about spatial position, intrinsic mean colour and magnitude, as 
well as pulsation properties of RR Lyrae stars. Using the stars’ 
coordinates, we removed variables caught at random phase in 


3. Synthetic Horizontal Branch modelling 

We computed synt hetic HB models wit h the code developed and 
fully described by jSalaris et al.| ( |20T^ toget her with the BaSTI 
library of scaled solar evolutionary tracks ( jPietrinferni et al. 
|2004| ). The use of scaled solar models with the same total metal¬ 
licity [M/H] of Carina stars is justified, since a-enhanced evo¬ 
lutionary tracks closely mimic scaled solar ones with the same 
tota l metallicity Z, in t he metallicity regime of this galaxy (see, 
i.e., |Salaris et al.|1993] ). 

We adapted the code, to allow for the large mass range on 
the Carina HB reaching -1.4 M©. Higher masses correspond to 
a population younger than 1 Gyr (at Carina’s typical metallicity) 
and can be ignored in our analysis, as their magnitudes in the 
helium burning phase are considerably brighter than the RC and 
are anyway a small contribution to the SFH of Carina. 

Briefly, for each each bin of the input SFH, a number of 
synthetic stars is generated spanning the whole range of age 
and metallicity inside the bin (employing a uniform probability 
within both age and metallicity bins), and with a value of [tr/Fe] 
given for that bin. For each star, the corresponding evolution¬ 
ary track is computed interpolating in mass and metallicity the 
tracks from the BaSTI grid. If the age of a star is greater than the 
age at the RGB tip, a specific amount of mass is removed and 
the position on the corresponding HB track is determined. 

Once the position of a star in the CMD is determined, its 
colour and magnitude are perturbed with a magnitude dependent 
Gaussian photometric error, as provided with BIO photometry. 
We then applied to the resultant synthetic CMD a reddening of 
E(B-V) = 0.06 ( jSchlegel et al.|1998] ) and a distance modulus of 
(m - M)o = 20.11, as used in dB14. Unless specified differently, 
we populated our synthetic CMDs with a considerably higher 
number of stars with respect to the observed CMD of Carina in 
order to minimize the Poisson error in the final model. 

The comparison of the model and the observed HB is a two- 
step process. First, an initial analysis is made by eye to see 
whether the colour extension and the main features of the HB 
are recovered. Then, after rescaling the total number of HB stars 
in the synthetic sample to the observed counterpart, we com¬ 
pare star counts and the mean colour and magnitude inside three 
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Fig. 4. Upper panel: Observed CMD of the HB region of the 
Carina dSph from BIO. Red crosses mark the position of RR 
Lyrae stars. The red boxes are those chosen for detailed compar¬ 
ison. The solid black line represent the ZAHB, extending up to 
a mass of 1.45Mo, for a metallicity of Z=0.0001. Lower panel: 
The synthetic HB for Carina. The red, blue and green dots mark 
stars with ages t > IGyr, AGyr <t< 1 Gyr and t < AGyr respec¬ 
tively. 

Table 1. The integrated RGB mass loss prescription 

as originally determined for the Sculptor dSph and used for the 
older Carina population, together with the one used for the inter¬ 
mediate age population in § |^ 


Metallicity range 

Sculptor 

Mo 

AMf°® 

t<lGyr 

Mo 

[M/P] < -1.8 

0.095 

0.048 

-1.8 < [M/P] < -1.3 

0.14 

0.07 

-1.3 < [M/P] 

0.16 

0.08 


boxes that encompass the RC (hereafter box RC), the red HB 
plus the RR Lyrae Instability Strip (IS - hereafter box R) and 
the blue HB respectively (hereafter box B -see Fig. [^. If the 
star counts are reproduced within one cr Poisson uncertainty 
and the difference in the mean photometric properties are within 
+ 0.01 mag, we consider the synthetic model to be a good match. 

We stress that we take into account only the uncertanties aris¬ 
ing from the Poisson distribution of the star counts. Another 
source of uncertainties is represented by the errorbar associ¬ 
ated with the star formation rate in every bin of our reference 
SFH. Unfortunately, including this uncertainty in our analysis is 
problematic, as the errors are not independent from each other. 
Indeed, the conservation of the total number of star and of the 
density distribution across the CMD required by the SFH fit¬ 
ting procedure introduces a correlation among the individuals 
errorbars, so a realistic modelling of the resulting uncertainty 
on the model HB star counts is unfeasible without the covari¬ 
ance matrix of the SFH solution that is not provided by dB14. 


Fig. 5. Upper panel: Observed (red) and synthetic (blue) star 
counts as a function of {B - V) colour in boxes B and R (left and 
right respectively). The bin size is 0.015 mag. Poisson errors on 
the observed star counts are also displayed. Lower panel: As the 
upper panel, but for the V magnitude. 


Undoubtedly, the issue of a proper inclusion of the SFH errors 
deserves to be addressed in future works. 


Clearly, theoretical HB models have also intrinsic uncertain¬ 
ties, that affect the predicted HB star counts as a function of 
colour and magnitude. The main source of uncertainty is re¬ 
lated tothetreatmentqft^ He-core convective mixing (see. 


e.g., |Cassisi & Salaris|[2013[ |Gabriel et al.||2Q14| [Spruit [2015 
and references therein) that affect evolutionary timescales, lu¬ 
minosities and morphology of the HB tracks in the CMD. The 
treatment of core mixing in the adopted BaSTI models includes 
semiconvection, as described in jPietrinferni et al.j ( |2004| ), and 
allows to reproduce the so-called R 2 -parameter, defined as the 
number ratio of asymptotic giant branch to HB stars, measured 
in a sample of galactic Globular clusters, that is very sen sitive to 
the treatme nt of core mixing during the HB phase (see jCassisi 
et al. 12003} for a thorough discussion) 


We also computed the pulsation period of synthetic stars in¬ 
side the IS and compared them with the observed period dis¬ 
tribution. We employ ed both the IS boundaries a nd the pulsa- 
tional equation from jPi Criscienzo et aT ( | 2004 | ), assuming a 
mixing length ml = 1.5 Hp. This choice reasonably matches the 
IS boundaries as inferred from the RR Lyrae colour distribution. 
The First Overtone periods {Ppo) are related to the Fundamental 
(P/ 7 ) by the relation log P /7 = log P/70 + 0.1 3 (|Di Criscienzo et al. 
[ 2004 ] ). 


As regards the stars in the so called OR zone, in which a RR 
Lyrae can pulse F, FO or double mode, we treated them as F 
pulsators. This is, in principle, a rough approximation but, given 
the small fraction of FO pulsators in the RR Lyrae sample (~ 
14%), we expect the number of FO pulsators in the OR zone to 
be very small. 
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Fig. 6. Observed (black) and synthetic (red) period distribution 
for RR Lyrae stars. The bin size is 0.05 day. Poisson errors on 
the observed counts are also displayed. 
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Fig. 7. Distribution of our HB synthetic stars in the age- 
metallicity plane. Upper panel: Synthetic generated employng 
our toy model SFH. Lower panel: Synthetic generated employ¬ 
ing the dB14 SFH. 


3.1. Results with the reference simulation 

As a star ting point, we empl oyed the mass loss law as deter¬ 
mined by [Salaris et al.| ( |2013| ) for the Sculptor dSph, that gives 
an increasing value for the integrated RGB mass loss with in¬ 
creasing metallicity, as described in Table Figure]^ shows the 
HB of Carina as well as the synthetic model computed from the 
SFH. Different colours mark different age ranges. For the sake 
of comparison, we show a synthetic HB with approximately the 
same number of stars as observed in Carina’s HB. 

We can notice two major differences between synthetic and 
observed HB. Firstly, what is manifestly poorly reproduced is the 
morphology of the RC, which is more disperse than the observed 
one, and with a “blue cloud” of stars that blend with the red end 
of the old HB. Some problems are not completely unexpected 
since Sculptor, for which the mass loss was determined to be de¬ 
pendent only on the metallicity, hosts only an old stellar popula¬ 
tion. Given the mass range of the stars on the RC, their lifetime 
on the RGB is considerably shorter with respect to the old pop¬ 
ulation. We can, therefore, naively expect the total mass loss to 
be different, plausibly lower for this population. We, thus, tried 
to modify the mass loss for the younger population, to reproduce 
the observations. However, in this way, the problem could only 
be mitigated but not totally solved. 

This can be explained as follows: the RC “blue cloud” above 
the HB is primarily composed of stars with Z < 0.0002. The 
black solid line in Fig. shows the Zero Age Horizontal Branch 
(ZAHB), extending to a mass of 1.45 M©, for Z = 0.0001, which 
is typical of the bulk of the metal poor population in Carina. The 
position on the ZAHB becomes redder as the mass increases un¬ 
til, over a certain mass threshold, the ZAHB turns toward bluer 
colours, without merging with the RC. Because of this, stars 
of that metallicity will always be bluer than the observed RC 
colour, regardless of the mass loss. 


Given that the synthetic HB is uniquely determined by the 
mass loss law and the input SFH, this problem is an indication 
that the properties of Carina stellar population recovered by the 
dB14 SFH model do not match exactly the true SFH. 

Considering now the HB within boxes B and R (the old HB, 
that comprises synthetic stars with age t > 7Gyr, displayed as 
red dots in Fig. |^, we can see that the total colour extension is 
nicely reproduced. This implies that, given the input SFH, the 
true mass loss for the old population cannot be drastically differ¬ 
ent from that inferred for the Sculptor dSph. Figure [^compares 
the colour and magnitude distributions of observed and synthetic 
HBs. As previously stated, the Sculptor-like mass loss law well 
reproduces the colour extension of the HB. Furthermore, the rel¬ 
ative star counts inside the two boxes are reproduced within one 
sigma. 

As for the mean colour and magnitude, the fit is not perfect. 
In box B, the computed mean colour differs by about 0.01 mag 
from the observed one, while the difference in magnitudes is of 
the order of 0.03 mag. It should be said that, in this region of the 
CMD star counts are low and a few stars can significantly alter 
the mean colour and magnitude. The synthetic distributions look 
very consistent within the observed error bars, suggesting that 
the difference in the means might be due to stochasticity. 

In box R the V magnitude mean value and the overall distri¬ 
bution are remarkably well reproduced. The difference between 
observed and synthetic mean V is of the order of the photomet¬ 
ric uncertainty. What are strikingly different are the star counts 
as a function of colour. The mean (B - V) colour of our synthetic 
stars differs by more than 0.05 mag from the observed one. This 
can be explained by looking at the colour distribution of the two 
populations. The observed star counts drop for colours bluer than 
0.4, and tend to increase again toward the blue HB. This “gap” 
around (B - V) = 0.3 is missing in the synthetic HB which is 
uniformly populated. 
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Fig. 8. Observed CMD of Carina (black dots) superimposed to 
the synthetic CMD computed from our toy SFH. Colour coding 
is the same as Fig.|^ 

Table 2. Age, metallicity and normalized Star Formation Rate 
of the toy SFH (see text for details). 



^min 

Gyr 

^max 

Gyr 

< [M/H] > 

SFR 

Burst 1 

3.0 

4.0 

-1.69 

0.12 

Burst 2 

4.0 

6.0 

-1.71 

0.39 

Burst 3 

8.0 

10.3 

-111 

0.29 

Burst 4 

10.5 

12.5 

-2.23 

0.2 


We highlight here, that the procedure followed by B10 to se¬ 
lect Carina member stars is less efficient for intermediate colour 
objects (0.4 < (5 - V) < 0.6) and a residual amount of field stars 
is left in the CMD, as it can be seen in Fig. 0it is interesting 
to note that this contamination affects the star counts in box R 
and may be partly responsible for some of the observed stars at 
the red side of the gap, with {B -V) > 0.4. A solid way to deal 
with this effect is hard to find, but we estimated that effect on 
the colour distribution to be of the order of 15-20% of the typi¬ 
cal star count. We conclude that, despite the contamination, the 
observed “gap” feature is true. 

The discrepancy in the distributions between model and data 
in the red box, which also includes the IS, can also be seen in 
Fig. which shows the period distributions for the RR Lyrae 
pulsators. The observed RR Lyrae population has a period dis¬ 
tribution strongly peaked around P ~ 0.6 J and is mainly com¬ 
posed of F pulsators, whereas the synthetic population has a 
much broader F period distribution and tends to strongly overes¬ 
timate the FO population. 

This can be understood by looking at the relation employed 
for computing the periods. The driver of the pulsational period 
is the effective temperature. The periods and the colour distri¬ 
bution are thus closely related. Both the lack of FO pulsators 
and of F short periods is explained by the observed “gap” in the 



V V V 


Fig. 9. Same plot as Fig. [^except for our toy SFH. The rightmost 
panel represents the RC box. The scale on the ordinate axis is 
set to 2, 2 and 10 counts per tickmark for the three upper panels 
(from left to right) and 1, 2 and 5 counts per tickmark for the 
lower panels. 

stellar distribution at the level of the IS. Our synthetic, more ho¬ 
mogeneously populated, HB naturally covers a broader range of 
periods and predicts too many stars in the FO zone of the IS. 

To overcome this problem we tried to change the mass loss 
law. An analysis of age and metallicity distributions along the 
HB revealed that our synthetic IS is populated by a broad range 
of ages and metal contents. Consequently, changing the amount 
of mass loss for a single value or a small range of metallicities 
(or ages) cannot create such a well-defined gap. On the other 
hand, changing the whole mass loss law will considerably affect 
the colour extension of the blue HB, which is extremely sensitive 
to small mass changes. The conclusion is that the observed gap 
in the HB cannot be reproduced by varying the mass loss law, 
except by adopting a very fine-tuned ad hoc dependence on both 
age and metallicity, which we see as physically hard to justify. 

We have therefore investigated the SFH as a possible cause 
of the discrepancies encountered, to assess the additional con¬ 
straints that can be made on the Carina dSph SFH using the HB 
modelling. 

4. Implications for the Star Formation History 

A complete characterization of the SFH based on the HB is 
not feasible, due to the uncertainty on the RGB mass loss. 
Nonetheless, as seen in the previous section, there are issues in 
our modelling, i.e. the presence in our synthetic CMD of the RC 
blue cloud and the missing HB gap around -V) ~ 0.3, that 
cannot be overcome by simply changing the mass loss. The only 
other option to solve this problem should therefore be connected 
to the SFH solution we are employing. 

We next investigated which component of the SFH could be 
responsible for these problems. The lower panel of Fig.|^shows 
the distribution of our synthetic HB stars in the age-metallicity 
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plane. The region of our synthetic HB wich correspond to the ob¬ 
served gap has been found to be mainly populated by two com¬ 
ponents which are enclosed within the two blue solid boxes: a 
group of metal poor stars with ages between 9 and 10 Gyr and a 

group of very old, more metal rich stars, with [MjH] -1.6. In 

particular, we note that the presence of a population of old metal 
rich stars, is also principally responsible for the broad synthetic 
RGB observed in Fig.[^ as these stars populate its reddest part. 

On the other hand, the RC blue cloud feature is caused, as 
noted in the previous section, by stars belonging to the interme¬ 
diate population and with [MjH] < -1.9 (red dashed box in the 
figure). 

This finding highlights how the HB morphology can point 
out “forbidden” ranges of age and metal content that would 
result in unobserved features. It also suggests that the input 
SFH has too wide a distribution in age and metallicity, whereas 
Carina’s true SFH is more confined to specific regions of the 
age-metallicity plane. To verify this hypothesis we built a “toy” 
bursty SFH. We tried to reproduce a synthetic HB as similar as 
possible to the observed one, trying at the same time to repro¬ 
duce qualitatively the RGB and the TO region. The best fit model 
is composed of four separate bursts whose properties are summa¬ 
rized in Table We employed a fiat probability distribution for 
the ages in the intervals given in the table, and a Gaussian [M/H] 
distribution with the listed average values and cr=0.1 dex. We 
adopted the Sculptor-like mass loss law for the old population. 
For the intermediate population, we assumed the RGB mass loss 
rate to be half that of the older population, to roughly take into 
account the shorter RGB lifetime (as an example, a O.SM© red 
giant, typical of a 12-13 Gyr old population, takes -1.4 Gyr at 
Z = 0.0001, and ~ 2.3 Gyr at Z = 0.001, to reach the RGB Tip 
from the TO. In contrast, the time taken by a IM©, typical of 5-6 
Gyr old populations, is - 0.7 Gyr and - 1.3 Gyr, respectively). 

Figure shows the resulting synthetic CMD superimposed 
on the observed one. Our toy model roughly reproduces the mor¬ 
phology of the TOs and the SGBs, as well as the RGB colour 
and the morphology of the HB. Figureshows the colour and 



Fig. 11. Left panel: Normalized MDF from Red Giants stars. 
The black histogram shows the MDF inferred from our toy 
model. The red shaded histogram represents the empirical dis¬ 
tribution determined by |Helmi et~ar] ( |2006| ) and employed by 
dB14. Central panel: as for the right panel but with the inclu¬ 
sion of intrinsic measurement errors in our MDF. Right panel: 
as for the central panel but also taking into account the observed 
spread in the [Ca/Fe] values of Carina stars (see text for details). 

magnitude distributions inside the B, R and RC boxes in Fig.|^ 
The relative numbers inside each box are reproduced within one 
sigma, and the mean values of colour and magnitude differ by 
at most 0.01 mag from the observed ones. Even the detailed dis¬ 
tribution looks consistent, although still with a few minor mis¬ 
matches . 

As an additional test, we looked at the RR Lyrae period dis¬ 
tributions, as shown in Fig.[^ The F period distribution is very 
peaked and looks similar to the observed one. The ratio between 
F and FO pulsators is qualitatively reproduced as well. 

It should be noted that this toy model is not intended to be the 
exact SFH solution, since we still lack strong constraints on the 
mass loss and we do not try to match the precise number density 
distribution across the whole CMD. Our purpose is to show how 
a more bursty SFH, in terms of age and metal content, is able 
to explain at the same time the morphology of the CMD and the 
detailed structure of the HB, and to suggest the age and duration 
of these bursts. 

The age-metallicity distribution of HB stars simulated with 
this model is shown in the upper panel of Fig.The comparison 
with the distribution in the lower panel derived from the dB14 
SFH allows us to identify differences and similarities between 
the two SFHs. The locations, in the age-metallicity plane, of all 
the star formation events of the toy model are roughly matched 
by those predicted by the dB14 solution, although the latter has 
a greater age spread, especially at old ages. Unsurprisingly, the 
major difference is the lack of metal poor stars for ages younger 
than 10 Gyr and of very old stars with [MjH] -1.6. 

In view of this comparison, we conclude that the discrepan¬ 
cies between the observed HB and the synthetic helium burning 
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population computed using the dB14 SFH are caused by two 
factors. First, the age resolution of the SFH naturally leads to 
a smooth HB, erasing any sharp substructure in the true SFH. 
dB14 test the age resolution of their method with a series of 
synthetic populations generated by a 10 Myr burst centered at 
different ages. Even in the best case, the recovered burst dura¬ 
tion would have been 500 Myr, that is the width of a bin element 
at old ages. It can be seen from their Fig.6 that in the old age 
regime, the recovered SFH is a Gaussian with a dispersion of 
the order of 1 Gyr. This issue is due to the strong degeneracy, 
in the age-metallicity space, of TO and RGB, coupled with the 
non-negligible photometric error at that magnitude level, for a 
system as distant as a dwarf galaxy. 

As a second point, the presence of spurious components in 
the dB14 SFH causes a fraction of the synthetic HB population 
to reside in regions of the CMD which are observed to be devoid 
of stars. The origin of these spurious components is not clear. A 
possible explanation is that any difference between the spectro¬ 
scopically measured MDF, which is used as a strong constraint 
in the SFH characterisation, and the true MDF may force the 
SFH determination algorithm to assign the wrong age to a frac¬ 
tion of stars in order to match the density distribution across the 
TO and the SGB regions. Regarding the Carina dSph, the prob¬ 
lem of metallicity is particularly thorny, with photometric and 
spectroscopic estimates, both high and mid resolution, not al¬ 
ways agreeing with each other (|Smecker-Hane et al. 


et al.l 2003} [Tolstoy et al.j|20Q3t [Bono et al.||2010t [Koch et af 
2006[ Lemake et al.|2012|. 


To investigate this point, we adopted the [tr/Fe] vs. [M/H] re¬ 
lation, recovered from the dB14 SFH, to compute the [Fe/H] val¬ 
ues of our model. We computed our MDF from RGB stars, down 
to 3 magnitudes below the RGB tip, which is the typical selec¬ 
tion criterion for DGs spectroscopic measurements. Inferring the 
metallicity distribution from the same region of the CMD is cru¬ 
cial, since, due to the varying evolutionary lifetime with metal¬ 
licity, stars end up on the upper RGB with a different MDF with 
respect to the original one on the MS. 

The left panel of Fig. pT] shows our synthetic MDF compared 
with the measurements from [Helmi et al.| ( |2006| ), which have 
been employed in the dB 14 analysis. In contrast with our sharply 
bimodal distribution, the measured MDF has a much broader 
distribution, which is unsurprising, given that the measurement 
errors naturally tend to smooth the underlying distribution. 

To check whether the two MDFs are consistent within the 
uncertainties, we convolved our MDF with a Gaussian error of 
0.1 de x, which is the typical uncertainty of the |Helmi et al.| 
( |2006|) m easurements. As can be seen in the middle panel of 
Fig. |11[ the bimodality of our MDF is still clearly noticeable 
and we conclude that the intrinsic measurements errors are not 
big enough to reproduce the observed MDF width which, in ad¬ 
dition, is fairly asymmetric. 

It is however important to notice that this medium reso¬ 
lution metallicity distribution has been inferred from the Call 
triplet (CaT). As [Battaglia et al.| ( |2008| ) showed, the equiva¬ 


lent width (EW) of this feature is sensitive not only to the iron 
abundance but also to the [Ca/H] value (as can be expected). 
Furthermore, the variation of other a-element abundances, espe¬ 
cially the [Mg/H] ratio, modifies the free electron density and, 
thus, the continuum opacity, indirectly affecting the measured 
EW. This effect has bee n already observed in globular clusters 
( jMucciarelli et al.|2012| ). 

The MDF adopted for the dB14 SFH has been evaluated us¬ 
ing the CaT-[Fe/H] calibration by jStarkenburg et alT] ( |2010| ). This 
relation has been determined assuming a fixed metal mixture 


with [Qf/Fe]=0.4 and [Ca/Fe]=0.25. As noted in the same work, 
at fixed iron abundance, a decrease of the [Ca/Fe] value of 0.25 
dex causes a smaller EW of the CaT feature, mimicking a value 
of [Fe/H] lower by ~ 0.2 dex 

Concerning Carina, several high resolution spectroscopic in¬ 
vestigations have found a large spread in the abundance of ele¬ 
ments like Ca and Mg, which has been explained as evidence of 
inhomogeneous mixing of the gas phase (jShetrone et al. 2003 [ 
[Lemasle et al.||201^ [Venn et ar||2012| [Fabrizio et al.||2015^ 
In particular, the bulk of the [Ca/Fe] measurements range from 
~ -hO.3 to slightly subsolar values. Using the CaT calibration 
without considering the variation in calcium abundance will then 
tend to un derestimate t he iron abundance. This has been already 
noticed by Venn et alT] ( |2012| ). 

A precise quantitative evaluation of this uncertainty is very 
difficult, since it requires detailed knowledge of several quanti¬ 
ties, like the trend of the mean [Ca/Fe] value and of its spread 
at varying metallicity, and the accurate EW dependence on the 
electron donor elements’ abundance. Nevertheless, as a first ap¬ 
proximation we took this effect into account by perturbing our 
MDF with an asymmetric uncertainty: we assumed that half of 
the stars had their [Fe/H] overestimated, with a dispersion of 
0.1 dex, and the other half had it underestimated, with a disper¬ 
sion of 0.25 dex. This roughly takes into account the Carina ob¬ 
served distribution of [Ca/Fe] around the fiducial [Ca/Fe] =0.25 
employed in jStarkenburg et al.] ( |2010| ) CaT-[Fe/H] calibration. 

The resulting distribution is shown in the right panel of 
Fig.[TT] The two distributions are in good agreement, suggesting 
that, at least qualitatively, an uncertainty of the order of 0.2-0.3 
dex toward lower metallicities is able to produce the broad MDF 
we observe. Assuming this to be the case, this issue could likely 
lead to spurious components in the best fit SFH, as the estimated 
high fraction of metal poor stars is included in the solution, with 
the algorithm accordingly distributing them in age to preserve 
the relative numbers between the two populations. 

5. Conclusions 

We have performed a detailed simulation of the HB of the Carina 
dSph using recent and detailed SFH estimates by dB14. We 
found that the overall colour extension of the old HB is well re¬ 
produced with integrated an RGB mass loss that ranges from 0.1 
to 0.14 Mq as a function of [M/H], in complete agreement with 
was has been found for the Sculptor dSph. This concordance 
hints that the mass loss law derived for the Sculptor dSph may 
reproduce the HB morphology in also other dSphs with similar 
metallicity. 

On the other hand we found some discrepancies with the de¬ 
tailed colour distribution along the old HB and RC that required 
a modification of the input SFH by dB14. We have then built a 
“toy” bursty SFH that reproduces well observed HB star counts 
in both V and {B - V), and also matches qualitatively the RGB 
and the TO region, as well as the pulsational properties of RR 
Lyrae stars. 

The comparison between the MDF of our bursty SFH and 
the one measured from the infrared CaT feature using the CaT- 
[Fe/H] calibration by jStarkenburg et al.j ( |2010| ) shows a quali¬ 
tative agreement, once taken into account the range of [Ca/Fe] 
abundances measured in a sample of Carina stars, that induces 
a bias of the derived [Fe/H] distribution toward too low values. 
Although we stress that Carina is one of the few Local Group 
galaxies where this issue can be so severe, extra caution should 
be taken when using the CaT to estimate [Fe/H] in system with 
poor or no constraints about the detailed metal mixture. 
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In conclusion, the work we have presented shows how the 
information contained within the HB can, in principle, be ex¬ 
tracted and interpreted to make predictions about the properties 
of the stellar population in DGs. The results on the Carina dSph 
illustrate that the inclusion of the detailed HB morphology in 
the CMD analysis is able to refine the SFH determination of re¬ 
solved stellar systems. Due to the strong sensitivity of its mor¬ 
phology and luminosity to the age and metallicity of the parent 
population, the HB is a powerful benchmark to rule out “forbid¬ 
den” ranges of age and metal content that would result in unob¬ 
served features on the CMD. This capability grants us additional 
resolution when pinpointing the SFH details of DGs, especially 
at early times. 

We aim to consolidate the synthetic HB modelling tech¬ 
nique, with the goal of making quantitative predictions on the 
properties of resolved stellar populations, consequently enhanc¬ 
ing the power of CMD analysis, especially in distant systems 
where the old MS-TO information is hard to acces. 
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